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FOREWORD 


Equilibrium  and  stability  properties  of  an  intense  electron  ring  located  at 
the  midplane  of  an  externally  applied  mirror  field  are  investigated  within  the 
framework  of  the  linearized  Vlasov-Maxwell  equations,  including  the  important 

influence  of  equilibrium  self  fields  and  an  applied  field  Bg^  in  the  toroidal 

direction.  It  is  assumed  that  the  ring  is  thin  and  that  v/y.  «  1,  where  v  is 
Budker's  parameter  and  y.mc^  is  the  characteristic  electron  Denergy.  Equilibrium 
and  stability  propertiesDare  calculated  for  the  choice  of  equilibrium  distribution 
function  in  which  all  electrons  have  the  same  value  of  energy  in  a  frame  of 
reference  rotating  with  angular  velocity  u)l  in  the  minor  cross  section  of  the 
ring,  and  a  Lorentzian  distribution  in  canonical  angular  momentum  PQ.  Negative- 
mass  and  resistive-wall  stability  properties  are  calculated,  and  a  ^closed 
dispersion  relation  is  obtained  for  the  case  where  the  ring  is  located  inside  a 
toroidal  conductor  with  finite  resistivity  and  minor  radius  ap  much  less  than 
the  major  radius  Rg.  One  of  the  most  important  features  of  the  stability  analysis 
is  that  the  negative-mass  instability  in  a  high-current  ring  can  be  stabilized 
by  equilibrium  self-field  effects  in  circumstances  where  the  self  fields  are 
sufficiently  intense.  Moreover,  a  modest  spread  A  in  canonical  angular  momentum 
can  stabilize  the  resistive-wall  instability. 


IRA  M.  BLATSTEIN 
By  direction 
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INTRODUCTION 


There  is  considerable  recent  interest  in  the  equilibrium  and  sta¬ 
bility  properties  of  intense  relativistic  electron  rings  with  applications 

1-4 

that  Include  high-current  betatron  accelerators.  In  a  conventional 

betatron  accelerator,  a  toroidal  electron  ring  is  confined  in  a  mirror 

or  betatron  magnetic  field.  Moreover,  the  beam  density  and  current  is 

limited  by  the  strength  of  the  betatron  magnetic  field.  In  the  modified 

1  2 

betatron  accelerator,  *  however,  an  additional  confining  magnetic  field 
Bqq  is  applied  in  the  toroidal  direction,  thereby  considerably  increasing 
the  limiting  beam  current  and  density.  In  the  present  article,  we  make 
use  of  the  Vlasov-Maxwell  equations  to  investigate  the  equlibrium  and  sta¬ 
bility  properties^  ^  of  an  intense  relativistic  electron  ring  confined  in 

1  2 

a  modified  betatron  field  configuration,  *  including  the  important  influ- 

pYf 

ence  of  equilibrium  self-field  effects  and  the  toroidal  magnetic  field  B  . 

U0 

The  analysis  is  carried  out  for  an  electron  ring  located  at  the  midplane  of 
an  externally  applied  mirror  field.  The  positive  ions  form  an  immobile 
(m^  -*■  «)  background  that  provides  partial  charge  neutralization.  In  addi¬ 
tion,  it  is  assumed  that  the  electron  ring  has  minor  dimensions  much  smaller 
than  the  major  radius  R_  (Fig.  1).  It  is  also  assumed  that  v/y.  <<  1, 

U  D 

2 

where  v  is  Budker’s  parameter  and  y  me  is  the  characteristic  electron 

D 

energy  associated  with  the  toroidal  motion. 

Equilibrium  and  stability  properties  are  calculated  for  the  specific 
choice  of  the  equilibrium  electron  distribution  function  [Eq.  (13)] 


«(H  -  u.  P.  -  YOC  ) 

_  p  y 

2’V  HP,  -  P0>2  +  a2] 
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where  H  is  Che  energy,  F  is  Che  canonical  angular  momentum  in  Che  azimuthal 

2 

direction,  P^  “  pp#  -  eS^p  /2c  is  the  canonical  angular  momentum  in  the 
small  cross  section  of  the  torus,  Rq  is  the  major  radius  of  the  ring,  and 
n^.  A,  Pq,  a>k  and  y  are  constants.  Although  H  and  PQ  are  exact  single-par- 
tical  invariants  for  the  equilibrium  configuration  illustrated  in  Fig.  1, 

Pf  is  an  approximate  invariant  whenever  the  toroidal  field  is  sufficiently 
strong  [Eq.  (8)]  and  the  electron  ring  has  circular  cross  section  (Appendix 
A)  with 


n  ■  —  and  a  ■  b  . 


Here  n  -  -[r(3/3r)£nB  (r,z)]/D  is  the  external  field  index.  Equili- 

uz 

8  9  u 

brium  properties  *  of  the  electron  ring  are  calculated  in  Section  II.  One 


of  the  most  important  features  of  the  equilibrium  analysis  is  that  the  max¬ 


imum  electron  density  trapped  in  the  ring  can  be  greatly  enhanced  by  the 
toroidal  magnetic  field  [Eqs.  (43)  and  (44)].  Moreover,  the  rotation  fre¬ 


quency  n>b  in  the  minor  cross  section  of  the  ring  plays  an  important  role  in 
determining  detailed  equilibrium  properties. 


The  formal  electromagnetic  stability  analysis  is  carried  out  in  Section 


III  within  the  framework  of  the  linearized  Vlasov-Maxwell  equations.  Neg¬ 
ative-mass  and  resistive-wall  stability  properties are  calculated 
for  elgenfrequency  u>  near  harmonics  of  the  cyclotron  frequency  wc2  ■  eB^  (R^ ,0) 

/y.mc  associated  with  the  axial  magnetic  field,  i.e.,  in  =  lw  .  The 
o  cz 

resulting  dispersion  relation  in  Eq.  (77)  is  obtained  analytically  for  a 


relativistic  electron  ring  located  inside  a  toroidal  conductor  with  finite 
resistivity  and  minor  radius  afi.  Equation  (77)  is  one  of  the  main  results 
of  this  paper  and  can  be  used  to  investigate  stability  properties  for  a 
broad  range  of  physical  parameters. 
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Negative-mass  stability  properties'"^ are  investigated  in  Section 
IV,  assuming  a  perfectly  conducting  wall  but  including  the  important  influ¬ 
ence  of  equilibrium  self  fields.  Introducing  the  parameter  [Eq.  (64)] 

2,2  2 
w  -  (*>cz/oig  -  i/yb  , 

we  find  that  [Eq.  (88)] 


is  a  necessary  and  sufficient  condition  for  instability.  Here,  g  is  a  geo- 

2 

metric  factor  with  g  *  (1  +  2ina  /a)  for  a  >>  a  and  g  •  4/A„  for  a  *  a, 

c  c  On  c 

where  J0(*0n)  "  °*  M0**0™**  -  cu^/2  +  0^/2)  [8*  “  U  ~  £)]  is  the 

radial  betatron  frequency-squared  for  a  circular  beam  with  a  *  b  and  n  ■  1/2, 

2,2 

u>pb  *  ^im5e  /t^®  re^at^v^-st^c  Plasma  frequency-squared,  and  f  is  the 

fractional  charge  neutralization.  Evidently,  from  Eqs.  (64)  and  (88), 
equilibrium  self  fields  can  have  a  large  influence  on  stability  behavior. 
Moreover,  the  negative-mass  instability  can  be  completely  stabilized  by  a 
sufficiently  large  spread  d  in  canonical  angular  momentum  P  . 


It  is  also  important  to  note  that  u  <  0  is  a  sufficient  condition  for  sta¬ 
bility  (Section  IV)  even  with  zero  spread  in  canonical  angular  momentum 
(A  ■  0).  The  condition  y  <  0  can  be  satisfied  provided  the  equilibrium 
self  fields  are  Sufficiently  strong  in  absolute  intensity.  For  example,  if 
f  *  0,  then  the  sufficient  condition  for  stability  (  y  <  0)  and  existence  of 

a  confined  equilibrium  [Eq.  (44)]  can  be  expressed  as  1  <  ^b^b^c  <  2  + 

2,-2 

U>  J  2  U  . 

CO  cz 

Although  a  perfect  conducting  wall  is  a  reasonable  assumption  in  many 
experiments,  we  expect  a  significant  modification  of  stability  behavior 
when  a  small  amount  of  resistivity  is  introduced  into  the  wall,  e.g.,  to 
stabilize  the  transverse  precession  in  the  minor  cross  section  of  the  ring. 
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In  this  regard,  in  Section  V,  we  investigate  the  resistive-wall  instability^ 


assuming  that  y  <  0  and  the  negative-mass  instability  is  absent.  It  is 
shown  that  the  resistive-wall  Instability  can  also  be  stabilized  by  a  modest 
spread  A  in  the  canonical  angular  momentum. 
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EQUILIBRIUM  CONFIGURATION  AND  BASIC  ASSUMPTIONS 


A.  Basic  Assumptions 


The  equilibrium  configuration  is  illustrated  in  Fig.  1.  It  consists 


of  a  relativistic  electron  ring  located  at  the  midplane  of  an  applied 
focussing  (mirror)  field  (r.z)^  +  Bg^Cr.z)^.  In  addition, 
the  electron  ring  is  located  inside  a  toroidal  conductor  with  minor 
radius  ac .  The  applied  toroidal  magnetic  field  Bqq  eQ ,  where 


‘  B 


ext 

08 


(1) 


together  with  the*  mirror  field  act  to  confine  the  ring  both  axially 
and  radially.  Here,  er,  e0,  and  are  unit  vectors  in  the  r-,  0-, 
and  z-directions ,  respectively.  The  equilibrium  radius  of  the  ring  is 
denoted  by  Rq,  and  the  minor  dimensions  of  the  ring  are  denoted 
by  2a  (radial  dimension)  and  2b  (axial  dimension).  In  addition  to  the 
cylindrical  polar  coordinates  (r,0,z),  we  also  introduce  the  toroidal 
polar  coordinate  system  (p,$,0)  illustrated  in  Fig.  1  and  defined  by 


r  -  Rq  =  r'  -  pcos$  , 

z  *  p  sin<t  , 


(2) 


where  p  is  measured  from  the  equilibrium  radius  Rq.  The  electrons 

composing  the  ring  undergo  a  large  orbit  gyration  in  the  external 

0 

mirror  field  with  characteristic  mean  azimuthal  velocity  VQ  ■  B^c  in 
the  positive  0-direction.  The  associated  ring  current,  which  is  in  the 
negative  0-direction,  produces  a  self-magnetic  field  $q(£)  that 
threads  the  ring  in  the  sense  indicated  in  Fig.  1.  This  self-magnetic 
field  acts  '  a  focuss-  \g  field  which  tends  to  confine  the  ring  electrons 
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both  axially  and  radially.  The  electron  ring  is  assumed  to  be 

partially  charge  neutralized  by  a  positive  ion  background.  The  excess 

electrons  form  a  potential  well  for  the  ions.  For  the  electrons, 

however,  the  self-electric  field  is  defocussing,  i.e.,  in  the  direction 

of  increasing  the  minor  dimensions  of  the  ring. 

To  make  the  theoretical  analysis  tractable,  we  make  the  following 

simplifying  assumptions  in  describing  the  electron  ring  equilibrium  by 

8  9 

the  steady-state  Vlasov-Maxwell  equations.  ’ 

(a)  The  positive  ions  form  an  immobile  (n^  ■+  •)  partially  neutralizing 
background.  The  equilibrium  ion  density  n^(r,z)  and  electron  density 
n^(r,z)  are  assumed  to  be  related  by 

n?(r,z)  -  fn?(r,z)  ,  (3) 

i  e 

where  f  *  const,  is  the  fractional  charge  neutralization. 

(b)  The  minor  dimensions  of  the  ring  are  much  smaller  than  its 
major  radius,  i.e.. 


a,b  <<  Rq  .  (4) 

To  further  simplify  the  analysis  we  also  assume  that  the  minor  cross 
section  of  the  ring  is  circular  with 

a  -  b  ,  (5) 

which  is  consistent  (strictly  speaking)  provided  the  external  field 
index  n  satisfies 

n  -  1/2  ,  (6) 

myf 

where  n  ■  -[rainB^  (r,z)/3r]  (Rq>o)  • 
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(d)  Consistent  with  Eq.  (4),  it  is  also  assumed  that  the 

transverse  (r,z)  kinetic  energy  of  an  electron  is  small  in  comparison 

2 

with  the  characteristic  azimuthal  energy  Y^mc  ,  i.e., 


lyTm  (Pr  + 


(?) 


where  p0  £  y^b^c  is  the  characteristic  azimuthal  momentum. 

(e)  The  maximum  spread  in  canonical  angular  momentum  6P„  is 
assumed  to  be  small  with  J  <SPQ |  <<  Yb®8bcRg  and 


5P„ 


Pu 

P0 


(8) 


2  2  1/2 

Here,  px  *  (pr  +  p^)  is  the  characteristic  transverse  (r,z) 
momentum,  p0  £  y^^S^c  is  the  characteristic  azimuthal  momentum, 
and  Bq  =  Bq^CRq.O)  and  B.,  -  Bq^CRq.O).  Moreover,  5Pg  =  PQ  -  PQ,  where 
P0  is  the  canonical  angular  momentum,  andPg  *  const,  is  the  average 
canonical  angular  momentum  of  the  electrons  composing  the  ring.  As 
shown  in  Appendix  A,  the  inequality  in  Eq.  (8)  (which  is  satisfied 
provided  the  spread  in  canonical  angular  momentum  £P  is  sufficiently 
small  and/or  the  toroidal  field  B  is  sufficiently  strong),  together 

•  O 

with  Eqs.  (5)  and  (6)  are  sufficient  to  assure  that  the  canonical 

*  2 

angular  momentum  P.  •  op.  -  (eB„/2c)o  defined  in  the  small  cross- 
9  9  9  - 

section  of  the  torus  is  a  good  approximate  invariant. 

(f)  It  is  further  assumed  that 


v 

Yv 


A-lLl  , 

2irR0  me2  <  ’ 


(9) 


2  2 

where  v  »  (Ng/2irRg)(e  /me  )  is  Budker’s  parameter,  N£  is  the  total 

2  2 

number  of  electrons  in  the  ring,  and  e  /me  is  the  classical  electron 
radius.  While  Eq.  (9)  assures  that  the  equilibrium  self  fields 


,  and  B*  are  weak  in  absolute  intensity  (in  comparison 
Or  u  z  o  z  or 


u 
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with  Bq  ‘  ),  depending  on  Che  beam  density  we  will  find  that  the 

self-field  gradients  can  be  large  and  have  a  correspondingly 

large  Influence  on  particle  orbits  in  the  equilibrium  field 
8  9 

configuration.  ’ 


B.  Self-Consistent  Vlasov  Equilibrium 

For  azimuthally  sysnnetric  equilibria  (3/39  -  0)  with  both  r  and 
z  dependence,  there  are  two  exact  single-particle  constants  of  the  motion. 
These  are  the  total  energy  B, 

H  -  (m2c4  +  c\2)V2  -  e*Q(r,z)  ,  (10) 

and  the  canonical  angular  momentum  P  , 

0 

♦ 

pe  "  r^pe  “  t  Ag(r*2>l  »  (ID 

where  £  *  ymv  is  the  mechanical  momentum,  $0(r,z)  is  the 
equilibrium  electrostatic  potential,  -e  is  the  electron  charge,  c  is 
the  speed  of  light  in  vacuo,  m  is  the  electron  rest  mass,  and  A?(r.z)  » 
Aqq  (r,z)  +  Ag0(r,z)  is  the  9-component  of  the  vector  potential 
for  the  total  (external  plus  self)  equilibrium  magnetic  field.  Without 
loss  of  generality,  we  assume  $q(Rq,0)  »  Aq0(Rq,O)  ■  0  in  Eqs.  (10) 
and  (11).  Within  the  context  of  Eqs.  (3)  -  (9),  it  is  shown  in 
Appendix  A  that  Che  canonical  angular  momentum 

P$  “  pp«  ”  2c  ®9P  ’ 

gxt 

in  the  plane  perpendicular  to  Che  toroidal  magnetic  field  B^0  is  an 
approximate  single-particle  invariant  for  a  thin  circular  beam  with  a  -  b 
and  n  •  1/2.  Here,  p  is  the  mechanical  momentum  in  the  ^-direction 
[Fig.  1],  and  p2  •  (r  -  Rq)2  +  z2  is  defined  in  Eq.  (2).  In  obtaining 
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Eq.  (12),  the  toroidal  magnetic  field  b' 


ext 

’06 


B0Rq/t  has  been 


approximated  as  uniform  over  the  minor  cross  section  of  the  ring. 

Any  distribution  function  f^  that  is  a  function  only  of  the  single¬ 
particle  constants  of  the  motion  in  the  equilibrium  field  configuration 
satisfies  the  steady-state  (3/3t  «  0)  Vlasov  equation.  For  present 
purposes,  we  consider  the  electron  distribution  function  specified  by 


0  5(H  "  “bP*  "  *mc  } 

fb(H’P*’V- 


2 

2ir  Yjjti 


2  2 
<*9  “  V  +  * 


(13) 


where  n^  »  o^(Rq,0)  is  the  electron  density  at  the  equilibrium 

radius  (r,z)  *  (Rq,0),  id^  *  const,  is  the  angular  velocity  of  mean 

rotation  in  the  ^-direction,  A  is  the  characteristic  spread  in  the  canonical 

angular  momentum  P  ,  and  y  is  a  constant. 

9  * 

The  equilibrium  Poisson  equation  can  be  expressed  as 

.2  \ 


i5rr»r+“l2jVr's> 


(14) 


4ue(l 


:i-f) 


• 


Furthermore,  the  6-component  of  the  V  *  B“(x)  Maxwell  equation  can  be 
expressed  as 


-  -  -  r +  ^|A;e(r,z) 


3r  r  3r 


3r 


(15) 


^  '»  'XW  > 


2  22  -1/2 

where  v  ■  (p  /m)(l  +  jp  /m  c  )  is  Che  azimuthal  electron  velocity. 
0  8  ^ 

Consistent  with  the  thin  ring  approximation  [Eq.  (4)]  is  the 
requirement  that  the  r-z  kinetic  energy  be  small  in  comparison  with 
the  effective  azimuthal  energy  [Eq.  (7)] 
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Y0(r,z)mc 


2 

me  +  c 


e  0. 

r  +  T  Vr-s) 


12  11/2 


(16) 


2  2  2- 

defined  for  p,  ■  p„  +  p„  ■  0  and  P  ■  Prt,  and  that  the  spread  in 
T  Z  0  U 

canonical  angular  momentum  be  small  with  |6P0|  ■  jPQ  -  Pg|  <<  ybm6jjcRg. 

2  2 

Taylor  expanding  for  small  values  of  p^  +  p„  and  6P  -  PQ  -  Po,  we 

r  z  o  c 

find  that  the  total  energy  H  defined  in  Eq.  (10)  can  be  approximated  by^ 


(17) 


2  2 

pr  +  pz  2 

H  ~  2~~r,z)m  +  Vr>z)mc  “  e*0<r>‘> 

6 

V?(r,z)  (6Pfl)2 

+  -8__  (5P  )  +  —-S - J  , 

2yg(r,z)mr 

where  the  mean  azimuthal  velocity  of  an  electron  fluid  element  V^(r,z) 
is  defined  by 


V^(r,z)  -  (/d3pv  f?)/ (/d3pf?) 


(18) 


a  [ (PQ/r)  +  (e/c)A0(r,z)]/[y0(r,z)m]  . 


For  future  reference,  we  define  the  characteristic  azimuthal  energy  Ybmci 
of  a  beam  electron  by 


2  2 
Ybmc  =  y0(Ro»O)®c  » 


(19) 


and  choose 


eB 


P«  - 


*  Jl 


2c  *0 


(20) 


A0(RO,°)  -  \  , 


where  =  b“C(R0,0). 
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The  mean  equilibrium  radius  Rq  of  the  ring  is  effectively 

8  9 

determined  from  the  condition  ’ 

p  17  Vr*z)  -  e  h  *o(r’z)]  (R0  0)  ‘  °-  (21) 

*  8  9 

Making  use  of  Eqs.  (16)  and  (21),  we  obtain  * 

ybm8^c2/R0  -  e[E°(R0,0)  +  0bB°(Ro,O)]  ,  (22) 

where  8b  =  Vg(RQ,0)/c,  E^(r,z)  -  -3<j>g/3r  is  the  radial  self  electric 

field  and  B^(r,z)  »  (1/r) (3/3r) (rA^)  is  the  axial  magnetic  field, 
z  0 

Equation  (22)  is  simply  a  statement  of  radial  force  balance  on  an 
electron  fluid  element  at  (r,z)  ■  (Rq,0).  Due  to  the  symmetry  of  the 
equilibrium  configuration  (Fig.  1),  it  is  clear  that  at  the  midplane 
(z  *  0),  the  axial  electric  field  and  the  radial  magnetic  field  are 
identically  zero.  Ue  therefore  note  that 
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we  Taylor  expand  Eq.  (17)  about  (r,z)  -  (RQ,0)  for  a  thin  ring, 


retaining  terns  to  quadratic  order  in  r'  »  r  -  Rn  and  z.  This  gives 
2 


7-10 


H  “  Y.nc  +  _ - 

b  2y.m 


2.2 

pr  pz  .  (.  r'  \ 

b-  cz  \  R0  /  I 


e 


1  2  2  2  2 

+  — +  2  v{v  ♦  «:*  >  . 

2Yb»Ro  * 

where  u^z  is  the  relativistic  electron  cyclotron  frequency  in  the 
axial  magnetic  field 


(27) 


eB 


cz  Yjj®c 


(28) 


and  Bz  =  Bq*  (Rq»Q).  Making  use  of  Eqs.  (16),  (25),  and  (26),  and  the 

expressions  for  the  equilibrium  potentials  given  in  Eqs.  (A. 6)  and 

2  2 

(A. 9),  we  find  that  u  and  u  can  be  expressed  as 

r  z 


■  “>?_(!  -  n)  +  TTT  [3v  “  (1-01  , 


cz 


pb  a+b  lPb 


and 


2  2  2  a  r„2  .  ... 

“z  "  “cz11  +  >  77b  tSb  ~  (1'f)]  ’ 


(29) 


(30) 


0  2  2 

where  Bb  *  Vg(R0’<^c’  “pb  *  ^Ybm  is  C^e  re^-ativist*c  plasma 

frequency-squared,  &b  •  n£(Rg,0)  is  the  beam  density  at  (r,z)  ■  (RQ,0), 


and 


_ext,  \ 
kB0z  (r’z) 


3  _ext,  . 
17  B0z  (r’z) 


(31) 


(R0,0) 


is  the  external  field  index. 

In  the  remainder  of  this  article,  the  equilibrium  and  stability 

analysis  is  restricted  to  circumstances  where  the  characteristic 

spread  &  in  canonical  angular  momentum  P  is  sufficiently  small  that 

8 
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A  «  Yb*wra  .  Ybt^raR0  . 

cz 


Within  the  context  of  Eq.  (32),  the  6 P  contributions  in  Eq.  (27) 

are  negligibly  small  in  comparison  with  the  final  term,  and  the 

total  energy  H  can  be  approximated  by 

2  2 

2  .  pr  pz  ,1  ,  2  ,2  .  2  2, 

~  +  -r  y.  m(tu_r  +  ut  z  )  . 


H  »  y.  me  + 


b“"  2y.  m  T  2  TbmvV  T  “V 


Referring  to  Eqs.  (12),  (13),  and  (33),  and  making  use  of  the  toroidal 

polar  coordinate  system  (p,$,8)  defined  in  Eq.  (2)  and  Fig.  1,  the 

combination  H  -  can  then  be  expressed  as 

p2 

H  ‘  “bp*  "  Ybmc2  +  TTZ  +  »  (34: 

u 

2  2  2 

where  p,  =  p  +  (p  -  Yv®UvP7  Is  the  transverse  momentum-squared 

p  $  b  b 

in  a  frame  of  reference  rotating  with  angular  velocity  iDfc  -  const, 
about  the  toroidal  axis,  the  envelope  function  4/(r',z)  is  defined  by 

>j)(r'  ,z)  -  j  Yb“(^yr'2  +  ft2 z” )  ,  (35] 

„2  -2 

where  r'  -  r  -  Ro,  and  fir  and  are  defined  by 

«2  2  2 
flr  "  “b“ce  '  “b  +  “r  ’ 


»2  2.2 

flz  "  wb“c9  "  “b  +  “z  * 

2  2 

where  w  and  w  are  defined  in  Eqs.  (29)  and  (30),  and 
r  z 


‘’ce  YvmC  ’ 


is  the  relativistic  cyclotron  frequency  in  the  applied  toroidal  magnetic 
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a  ext 

field  B0  -  Bq0  (Rq,0).  Substituting  Eqs.  (34)  and  (35)  into  Eq.  (13) 
and  evaluating  the  electron  density  profile  n^(r,z)  »  Jd^pf^, 
we  find 


nj(r,z)  -  t^U 


(■•(?•«] 


(38) 


where  a  »  [2(y  -  yfe)c  /rbSr3  and  b  ■  [2(y  -  Y^c2^!!2]1  *  0(*) 

is  the  Heaviside  step  function  defined  by 

1  ,  x  >  0  , 

0  ,  x  <  0  . 


U(x) 


Note  from  Eq.  (38)  that  the  electron  density  is  constant  (nb)  inside 

2  2  2  2 

the  beam  cross  section  r'  /a  +  z  /b  <  1,  and  identically  zero  outside. 

As  indicated  earlier  in  this  section  and  in  Appendix  A,  P$  is  a 

good  invariant  (strictly  speaking)  provided  Eq.  (8)  is  satisfied  and 

provided  the  beam  is  circular  with  a  »  b  and  n  *  1/2.  From  Eqs. 

(29),  (30),  and  (36),  we  therefore  find  w2  *  w2  *  ui2  and  02  ■  02  *  fi2 

r  z  p  r  z  p 

where 


and 


1  2 

T  w 
2  cz 


«; >1 


-  (i-f))  , 


«;■ 


“b“c8 


%  + 


(39) 

(40) 


Moreover,  the  minor  radius  of  the  ring  is  given  by 

a  -  b  -  [2  (y  -  Vc2/YbnB]1/2  *  (41) 

2 

For  exis cence  of  the  equilibrium ,  both  y  >  y,  and  ft  >  0  are  required. 

D  p 

2 

The  condition  0  >  0  can  be  expressed  in  the  equivalent  form  [Eqs.  (39) 

P 

and  (40)] 


“b  <  "b  4  “b  ’ 


(42) 
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where  Che  frequencies  u~  are  defined  by 


1  t 


2o» 


1  + 


cz 


u 


,  '1/2’ 


(43) 


c0  c9 

The  requirement  that  the  radical  in  Eq.  (43)  be  real  determines  the 


maximum  allowable  equilibrium  beam  density  for  given  values  of  f,  B? , 

D 


and  “ce- 


For  example,  if  f  ■  0,  then 


(44) 


2  2 

is  required  for  existence  of  the  equilibrium.  For  >>  w^, 

Eq.  (44)  is  identical  to  the  result  obtained  by  Sprangle  and 
Kapetanakos^  within  the  context  of  a  model  that  examines  the  stability 
of  single-particle  orbits. 

Finally,  we  conclude  this  section  by  noting  that  Eqs.  (13),  (32), 

(34),  and  (35)  can  be  used  to  evaluate  a  variety  of  equilibrium 
8  9 

properties  *  for  a  circular  beam  with  a  *  b  and  n  *  1/2.  For  example, 
the  mean  rotational  velocity  -  (/d^pv^f^) / (/d^pf^)  of  an  electron 
fluid  element  about  the  toroidal  axis  is  given  by 


V 


(45) 


where  the  rotation  frequency  is  restricted  to  the  range  in  Eq.  (42). 
Moreover,  it  can  be  shown  that  the  equilibrium  pressure  tensor  in  the 
(p,$)  plane  perpendicular  to  the  0-direction  is  isotropic  with  the 
perpendicular  pressure  P$(p)  *  n^(p)l2(p)  given  by 


n°(P)T°(p)  -  2*  L 


.  p0  +  <p»  -  .0 
dpe  - - £ 


b  * 
(46) 
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where  TA(p)  is  the  effective  transverse  temperature  profile. 


s.  .  1  .2  2  1  2  2 

T*  =  I  VV  ■  2  ^““ce'L  * 


and  substituting  Eq.  (13)  into  Eq.  (46)  gives 


T?(p)  -  Tx(l  -  p2/a2)  , 


Defining 

(47) 

(48) 


for  0  <  p  <  a.  In  Eq.  (47),  r^  is  the  characteristic  thermal 


Larmor  radius  of  the  ring  electrons  in  the  azimuthal  magnetic  field  B  . 

6 

Making  use  of  Eq.  (47)  to  eliminate  ft.  in  Eq.  (40)  in  favor  of  rT , 

P  L 


we  solve  Eq.  (40)  for  the  rotation  frequency  and  obtain 


.±  .  C0 

“b  “  “b  =~<1± 


2iu 


2<u 


2r, 
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1  + 


-s“-T*a 


C0 


ce 


1/2' 


(49) 


which  relates  u>b  to  the  Larmor  radius  rL»  The  two  signs  (±)  in  Eq.  (49) 
represent  fast  (+)  and  slow  (-)  rotational  equilibria.  In  order  for 
the  equilibrium  to  exist,  the  ratio  2rL/a  in  Eq.  (49)  is  restricted 
to  the  range 


<  1  - 


,  2u> 

(1  .  f  .  s2)  +  _« 
“ce 


(50) 
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LINEARIZED  VLASOV-MAXWELL  EQUATION'S 

In  this  section,  we  make  use  of  the  linearized  Vlasov-Maxwell 
equations  to  investigate  stability  properties  of  the 
equilibrium  ring  configuration  discussed  in  Sec.  II.  In  the  stability 
analysis,  a  normal-mode  approach  is  adopted  in  which  all  perturbed 
quantities  are  assumed  to  vary  according  to 


6\ii(g,t)  -  6$ (^c)exp (-iwt) 

-  $Mp)exp[i(£6  -  (lit) J  , 


2  2  1/2 

where  Imu  >  0,  p  i  [(r  -  R^)  +  z  ]  ,  and  3/3$  *  0  is  assumed 

(Fig.  1).  Here,  a>  is  the  complex  oscillation  frequency  and  l  is  the 
toroidal  harmonic  number.  Integrating  from  t'  ■  -®  to  t'  *  t  and 
neglecting  initial  perturbations,  we  find  that  the  perturbed 
distribution  function  can  be  expressed  as  6fjj(x,^,t)  *  5ffa  (^,^)exp(-iwt) , 
where 


r° 

6fb(x,jg)  «  e  dr  exp  (-iwt) 


+  £*'  *  «£<*’>} 


(51) 


3 _  e0 ,  ,  i  •> 

3^'  fb(£  ’■? 


In  Eq.  (51),  x  •  t'  -  t,  5£(£)  and  6£(£)  are  the  perturbed  electro¬ 
magnetic  field  amplitudes,  and  the  particle  trajectories  ^'(t')  and  v'(t') 
satisfy  d^c'/dt'  ■  }/'  and  djj'/dt'  «  -e[£^  +  v'  x  £^/c]  with  "initial" 
conditions  x' (t'  -  t)  ■  x  and  v'(t'  -  t)  -  v.  The  Maxwell  equations 

'v  %  'V 

for  6E(ic)  and  6BQc)  are  given  by 


7  x  £E 


iui 

c 


7  X  5B 


Uir 


(52) 

■ 

(53) 

.  t  ■  **'  i,  £  ./**** 
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where  the  perturbed  curr-nt  5 J  is  related  to  and  dE  by  -Ik:  rk 


(4ir/c)6Jp  -  (iu/c)6Ep.  We  also  assume  that  the  toroidal  conductor 
(radius  P  *  a£)  has  large  aspect  ratio  with 

ac  <<  V 


(60) 


Taking  the  8-component  of  Eq.  (57),  we  obtain  the  eigenvalue  equation 
for  6E0, 


»t?clr  +  (vk2)  Sk6  •  4’(lkss  -  7  sie)  • 


(61) 


2  2  2  2  -1 

where  k  *  Z/Rg,  and  we  have  approximated  V  -  -  Z  / r  *  p  (3/3p)  x 

2 

(p3/3p)  -  k  in  Eq.  (61).  It  is  further  assumed  that  Rew  -  satisfies 

JJ  a.  Zu>  ,  and  that  the  waves  are  far  removed  from  resonance  with  the 
r  a.  cz’ 

transverse  (r,z)  motion,  with 

+  2 
“b 


id)  -  Z<o 


cz 


)- 


(62) 


where  <n£  are  the  characteristic  (r',z')  orbit  oscillation  frequencies 
about  (r,z)  -  (Rq,0)  (Eq.  (43)  and  Appendix  B],  To  lowest  order, 
consistent  with  Eq.  (62),  it  is  shown  in  Appendix  B  that  the  azimuthal 
orbit  is  described  by  [Eq.  (B.14)] 


0' 


0  +  <«„  -  u6P0/YboRo)T  , 


(63) 


where  t  -  t'  -  t. 


to  , 

cz  1 

2  "  2  ’ 

“B  Yb 


(64) 


2  2  2  2 

and  (o  -oi  /2  +  (to  ,/2)[8.  -  ( 1—  f )  3  for  a  circular  beam  with  a  -  b  and 
g  CZ  po  D 

n  •  1/2. 

From  Eqs.  (51),  (58),  and  (59),  we  note  that 
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«%ce'> 


«*  * 


exp 


•*»  <it9'>  +  f  (vf?iVc')  +  T  9Vp,))]  5e  <95> 

t 

(!«•)[«(»•)  ♦^r*v,,',)]j>  • 


Estimating  the  size  of  terms  proportional  to  <$E  and  35E./3p*  in  Eq.  (65), 

P  6 

and  making  use  of  the  assumption  that  the  wave  perturbations  are  far  re¬ 
moved  from  resonance  with  the  transverse  (r’»  z*)  motion  [Eq.  (62)  and 
Appendix  B],  it  is  straightforward  to  show  that  the  terms  proportional  to 

A  ^  A  A 

36E  /3p*  and  6E  in  Eq.  (65)  can  be  neglected  in  comparison  with  6Ene. . 

9  P  4  9^0 

Moreover,  making  use  of  3f^(H,  P.,  P0)/3p.  ■  v„3f^/3H  +  r3f®/3P. ,  we  find 

D  V000D  D  0 

A 

that  5f,  (x,  p)  in  Eq.  (51)  can  be  approximated  by 

D  i\. 


3f, 


4fb(x,]^)  -  exp  (.±ld)f  dT(RQ  +  p'cos$’)6E9(p ')  exp  [ii(9T  -  9)  -  iuix] 

0  OJ 


(66) 


+  e 


3f“  r  0 

■jg  exp  (±l&y  dtVgdEg (p * )  exp  [ii(9’  -  9)  -  iut]. 


where  0'(t)  is  defined  in  Eq.  (63),  and  use  has  been  made  of  the 

fact  that  3f°/3P0  and  3f°/3H  are  independent  of  t’.  In  Eq.  (66), 

it  follows  from  Eq.  (63)  that  v'  -  (R  +  p'cos$')(w  -  ydP./y.mR^). 

0  u  cz  0  b  0 

To  simplify  Eq.  (66),  we  approximate  R^  +  p'cos$’  =  RQ  to  lowest 
order,  and  Taylor  expand  SE^p")  -  6E0(p)  +  [36E0  (p)/3p  ]  (p  •  -  p)  +  ... 
locally  about  p’  *  p.  For  present  purposes,  we  also  approximate 

•  A 

6E0(p')  *  6E0(p),  which  is  a  good  approximation  provided  r^  <<  a 
[Eqs.  (49)  and  (50)].  Carrying  out  the  t'  integration  in  Eq.  (66) 
then  gives 
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6f,(x,p)  =  ieexD  (  U'-  )■>  K.  (o  )R_. 
b  *V  0  u 

•af°  2  3f{\ 

JT  +  (ucz  -  u5VYbmV  wt 

8  _ 

2 

(u>  -  iw  +  i.u6Pa/y,  mRn) 
cz  o  d  u 


(6 


'  ! 


To  complete  the  description,  we  evaluate  4ir(ik5p  -  iw^Jg/c2)  * 

-4ireikJ  d  pSf.  (1  -  up  /c  kym)  on  the  right-hand  side  of  Eq.  (61).  To 


the  required  accuracy,  (1  -  up  /c  kym)  can  be  approximated  in  the  inte- 
grand  by  (1  -  w8b/ck)  -  w5PQ/c  kybmR0,  where  tu  a  *w.z  -  kw^Rg  *  kBbc  is 
assumed.  Making  use  of  Eqs.  (56),  (61)  and  (67)  then  gives  the  eigenvalue 
equation 


p 

/  2  „\1 

* 

”  / 

\  SP„  1 

1 

3  3  .fu  .2) 

6EQ(p)  =>  4re2i5Ee(p)y*d3p 

w  \  u  8 

p 

3p  P  3p  +  ^2  -  k  j 

[v 

Tj  ck  J  ck  Yj^mR^c 

3f 


iP„ 


b+l  _y_l 

3P  \  cz  _2 

8  \  TvniR, 


(68) 


in  -  +  !lii<5P  /y  mR 

CZ  8  D  U 


Paralleling  the  analysis  in  Ref.  10,  the  integral  contribution 
proportional  to  jd^ p  ...  3fb/3H  gives  surface-charge  contributions  on 
the  right-hand  side  of  Eq.  (68)  proportional  lo  3nb(p)/3o  *  -n^^p-a) 

[Eq.  (38)].  Although  this  surface-charge  contribution  can  be 
retained  in  a  self-consistent  manner,  for  present  purposes  we  neglect 
the  term  proportional  to  3fb/3H  in  Eq.  (68)  and  retain  only  the  negative- 

mass  and  resistive  wall  effects  associated  with  3f?/3Pa.  This  is  a 

Si  o  8 

0 

valid  approximation  provided  the  beam  density  satisfies  Eq.  (44). 
Integrating  over  3fb/3P0,  Eq.  (68)  then  reduces  to 


J 


NSWC  TR  81*389 


ii  3_  +  /  4! - 

[p  3p  p  3p  \c2  j 


«E0(p) 


Xh0*>) 

-2j-  U  (a  -  p)5E0(p)  , 


(69) 


where  the  susceptibility  xb(<*>)  is  defined  by 


2,22 


xb<«> 


(w  - 


Vk  * _  U  .  jA  .  M(U  I  *"cz>1  , 

c  +  l|ukA|/Yb»R0)2  L  V  bck/  C2k2 


for  the  choice  of  equilibrium  distribution  function  f^  in  Eq.  (13).  Since 

«  *  ■  kBbc  is  assumed,  the  square-bracket  factor  in  the  preceding  ex- 

2  2 

pression  for  xv(“)  can  be  approximated  by  u(l  -  S.  )  *  u/y.  and  the  suscep- 

D  D  D 

tibility  reduces  to  the  approximate  expression 


-  u 


2  .2*2.  2 

“Pbk  a  \ 


(w  -  !m>cz  +  i|ukA|/YKmRn) 


b  0 


(70) 


For  the  choice  of  f°  in  Eq.  (13),  we  note  from  Eq.  (38)  that  the  density 
profile  n^(p)  ■  i^lKa  ”p)  is  constant  in  the  beam  interior.  Here,  U(a  -  p)  - 
+1  for  p  <  a  and  U(a  -  p)  ■  0  for  p  >  a,  so  that  the  contribution  on  the 
right-hand  side  of  Eq.  (69)  corresponds  to  a  body-charge  perturbation. 

Of  course,  the  eigenvalue  equation  (69)  must  be  solved  subject  to  the 
appropriate  boundary  conditions.  At  the  beam  boundary  p  ■  a,  these  condi- 

A  A 

tlons  are  that  6E0(p)  and  36E0(p)/3p  be  continuous.  At  the  conducting  wall 
(p  «  a£  >  a),  we  assume  that  the  wall  has  finite  conductivity  a.  The  boun¬ 
dary  condition  at  the  wall  can  then  be  approximated  by^ 


sV*c>  •  -  H  a  -  «  sV»c> 


|b^u  +  1) 


[£  S‘Ee(0)]  «, 


(71) 


0 


iw 
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In  obtaining  H  .  (  .  .  *•  i<  fcuen  mad-  of  Eq.  (59)  on  the  vacuum  side 

*  ,  2  2  0 
(where  6j^  =  0)  of  the  conducting  wall,  and  we  have  approximated  (  u  /?  k”)> 

(1  -  w^/c^k^)  ^  =  SfY?  for  u  t  l*i  “  k3,  c.  In  Eq .  (71) 
d  D  cz  b 


6  *  c/ (2irwcr) 


(72) 


is  the  skin  depth  in  the  wall. 

Equation  (69)  can  be  further  simplified  for  jx^l  £  1  and  w  =  ^Qz  m 

ku  Rn.  In  this  case, 
cz  0 


which  is  consistent  with  Eq,  (62).  Equation  (69)  can  then  be  approximated  by 


I  1_  „L  +  t2 
(p  3p  P  3p 


6E0(P)  -  0  , 


and 


where 


£fc*!7*v*)  -  °  • 


I2a2  ■  xb<«> 


0  <  p  <  a  , 


a  <  p  <  a 
-  c 


(73) 

(74) 

(75) 


and  Xj,(“)  is  defined  in  Eq.  (70).  The  physically  acceptable  solutions 
to  Eqs.  (73)  and  (74)  with  5E  (p)  and  36E  /3p  continuous  at  p  *  a 
are  given  by 


«,<»>  ■  ( 

1  AJ0(Tp) 

1 

,  0  <  p  <  a 

f  J0(Ta)  » 

^  ^(Ta) 

1  +  la  .  ,  In  ^ 

l  VTa)  *> 

(76) 


a  <  p  <  a 
-  c 


where  Jq(x)  is  c^e  Bessel  function  of  the  first  kind  of  order  zero, 
and  Jq(x)  *  dJQ(x)/dx.  Enforcing  Eq.  (76)  at  the  conducting  wall 
g  -  a£  gives  the  required  dispersion  relation 
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VT*>  -  Jj^(Ti)  ^  -  I  (1+D  ^  8jf  j]  -  0  .  (77) 

where  use  has  been  made  of  J^(x)  »  -xjg(x). 

Equation  (77)  Is  a  transcendental  equation  for  the  conplex 

eigenf requency  u  and  must  be  solved  numerically  In  the  general  case. 

There  are  two  limiting  regimes,  however,  where  Eq.  (77)  can  be 

simplified  analytically. 

2  2 

(a)  ac  »  a  or  >>  a£:  There  is  a  large  temptation  to  Taylor 


expand  Eq.  (77)  for 


I  T|  2a2  «  1. 


To  lowest  order  this  gives  the  approximate  dispersion  relation 

From  Eq.  (79),  it  is  evident  that  Eq.  (78)  is  a  valid  approximation 


only  if 


j  (i  +  ii"  r )-  r  C1+1)8b1,b  >>x 


is  satisfied.  Equation  (80)  typically  requires  a  large  conducting  wall 
radius  ac  »  a  ,  or  bJyJ*  »  a£  in  order  for  che  inequality  to  be  satis¬ 
fied.  Making  use  of  Eqs.  (70)  and  (75),  the  dispersion  relation  in  Eq. 

(79)  reduces  to 


(w  -  i«cz  +  i|wkA|/Yb«Ro)' 


-  ?  ft1  +  2‘“  7^)  '  I"  <1+i)  ebVb]- 

y,  L  c 
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i  '22 

Equation  (77)  also  supports  solutions  with  ] T j  a  >1  when  the 
inequality  in  Eq.  (80)  is  satisfied.  Making  use  of  Eq.  (80),  we  find 
that  the  dispersion  relation  in  Eq.  (77)  can  be  approximated  by 

J^Ta)  -  0  ,  (82) 

which  gives  T^a^  »  X?  ,  n  =  1,2,...  where  J.(X,  )  *  0  and  X,  is  the 
°  In  1  In  in 

n'th  zero  of  J^(x)  *  0.  Making  use  of  Eqs.  (70),  (75),  and  (82), 
we  obtain  the  dispersion  relation 


(<D  -  iu>cz  +  i  I  pkA  I  /YbmRQ) 


”  -liU 


,22.2 
2  k  a  \  . 
pb  .2 
in 


1,2,. 


(83) 


where  ®  • 

Equation  (83)  is  really  an  extension  of  the  dispersion  relation 
2  2  2  2 

(81)  from  the  regime  [ T J  Za  <<  1  to  the  regime  |X|  a  >1.  What  is  most 

important  to  note  is  that  Eq.  (83)  is  independent  of  wall  resistivity, 

whereas  Eq.  (81)  is  not.  Moreover,  the  characteristic  growth  rate  «i  * 

law  obtained  from  Eq.  (81)  is  larger  than  that  obtained  from  Eq.  (83). 

2 

Indeed  the  mode  is  rapidly  stabilized  for  increasing  X  . 

2  2 

(b)  a£  a  a  and  B^o  <<  In  circumstances  where  BbYb<S  «  a£  and 

the  beam  radius  is  approximately  equal  to  the  conducting  wall  radius  a  a 
the  coefficient  of  J^Ta)  in  Eq.  (77)  is  algebraically  small,  and  the  dis 
persion  relation  (77)  can  be  approximated  by 


J0(Ta)  -  0  ,  (84) 

2  2  2  i 

which  has  solutions  T  a  *  Xqb,  n  *  1,2,...,  where  XQn  is  the  n  th 

zero  of  Jn(\+  )  »  0.  Making  use  of  Eqs.  (70),  (75),  and  (84),  the 
0  On 

dispersion  relation  (84)  gives 
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(w  -  + 


ijukij  /yboR0)2 


-yw 


,2  2 

k  a 


pbv2X2 

V 


1,2 


On 


(85) 


where  ^o^On^  "  ®*  most  unstable  solution  associated  with  Eq. 

(85)  Is  the  fundamental  mode  with  n  »  1.  As  in  Eq.  (83),  the  growth 

2 

rate  Xmu  decreases  rapidly  for  increasing  Xqq.  Comparing  Eqs.  (81) 
and  (85),  we  note  that  placing  the  conducting  wall  radius  (ac)  near 
the  beam  radius  (a)  significantly  reduces  the  growth  rate. 


30 


NSWC  TR  81-389 


NEGATIVE-MASS  INSTABILITY 

Vie  now  make  use  of  the  dispersion  relations  (81)  and  (85) 
derived  in  Sec.  Ill  to  determine  negative-mass  stability  properties 
for  the  case  where  the  wall  is  perfectly  conducting  with  6=0. 
Introducing  the  geometric  factor 


If* 


2tn 


20- 


for  a  >>  a  , 
c 


(86) 


4 

2 

^On 


for  a  =  a  , 
c 


2  22  2  222  22 
and  approximating  k  -  id  /  c  *  k  (1  -  w^Rq/c  )  ■  k  /y^  for  w  = 

the  dispersion  relations  (81)  and  (85)  can  be  expressed  as 


(id  -  iu>cz  +  ijpkAj /y^ihRq)' 


(87) 


-poj 


2  iV« 

pb  2  4 

Yb 


Solving  Eq.  (87)  for  w,  we  find  that  the  necessary  and  sufficient 
condition  for  instability  (Im^  >  0)  is  that 


0  <  p  < 


JL.  f2f 

Yb  \  A 


(88) 


2  2  2 

where  v  *  (fl^ira  )(e  /me  )  is  Budker's  parameter  for  the  beam,  and 
g  is  defined  in  Eq.  (86).  When  the  inequalities  in  Eq.  (88)  are 


satisfied,  the  real  frequency 


Reu  and  growth  rate  £1 


Imu 


are  given  by 


lu> 


cz 


n 


i 


1/2  „ 

i£_f  (  k£5 L )  -  vA 

YbRoi'Yb  mR0c- 


(89) 


where  k  ■  i/Rq. 
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The  stability  results  in  Eqs.  (87)  -  (89)  can  be  investigated 
in  various  regimes  of  experimental  interest.  It  is  interesting  to 
note  from  Eq.  (87)  that  a  sufficient  condition  for  stability  (Inw  <  0)  is 
u>2  i 

V  -  -p  -  If  <  0  ’  (90) 

“8  Yb 

2  2  2  2 

where  <o&  «  «*»cz/2  +  <o»pb/2>  [6b  -  (l-f)l.  For  f  -  0,  the  condition  y  <  0 
[Eq.  (90)]  and  condition  for  existence  of  confined  equilibria  [Eq.  (44) ] 
can  be  expressed  as 


1  < 


2  2 

Y.  oi 

b  cz 


<  1  + 


2u 


£0 

2 

cz 


(91) 


•Evidently,  for  f  *  0,  the  inequality  in  Eq.  (91)  can  be  satisfied  pro- 
2  2 

ided  0)  /u>  is  sufficiently  large.  That  is,  the  negative-mass  instabi- 

pD  CZ 

lity  can  be  completely  stabilized  for  4  ■  0  provided  equilibrium  self¬ 
field  effects  are  sufficiently  strong. 

The  regime  where  y  >  0  is  perhaps  of  more  practical  interest. 

In  this  case,  making  use  of  Eqs.  (87)  and  (88),  the  instability 
is  completely  stabilized  whenever  the  inequality 


1  4  ' 

l\v_l 

'  8  \ 

Vrbm8bcR0/ 

1  T„l 

(92) 

is  satisfied.  It  is  evident  from  Eqs.  (86)  and  (92)  that  the  instability 

is  most  difficult  to  stabilize  when  a^  >>  a  and  the  geometric  factor  g 

is  large.  Even  in  this  case,  however,  only  a  modest  spread  4  in 

canonical  angular  momentum  is  required  for  stabilization.  As  a 

2  2  2  2 

numerical  example,  consider  the  case  where  u>  , /u>  <<  1  and  y  =  1  -  1/v,  *  6,  . 

pb  cz  d  b 

For  v/y^  *  1/20,  g  ■  5,  Yjj  “  5,  and  8b  =  1,  Eq.  (92)  predicts  stability 
for  (A/Ybm3bcRg)  >  1/10. 
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RESISTIVE  WAIT.  INSTABILITY 


In  this  section,  we  consider  the  dispersion  relation  (77)  and  its 
limiting  versions  [e.g.,  Eq.  (81)]  in  circumstances  where  vail  resistivity 
effects  (as  measured  by  6/ac)  play  an  important  role.  In  this  regard, 
since  a  spread  A  in  canonical  angular  momentum  has  a  stabilizing 
influence,  we  first  consider  the  most  unstable  case  with 


and 


0  , 


_2  2 
T  a 


■2,22 

(d  .  k  a 

V  pb  _ 

2  2 
y  (to  -  ito  ) 
b  cz 


(93) 


(94) 


We  further  consider  circumstances  where 


u  <  0  , 


(95) 


and  therefore  the  negative-mass  instability  is  absent 

2  2 

Assuming  weak  dissipation  in  Eq.  (77)  with  S^y^6/ac<< 
lation  frequency  *  Re<o  is  determined  from 

J0(TO  -  ■>!«,)  0  , 

and  the  growth  rate  *  Imio  is  determined  from 


(Sec.  IV). 

1,  the  real  oscil- 


(96) 


fii 


c 


JJ(Ia)  -  ir&l) 


(97) 


3ft 


Ta 


where  Ta  in  Eq.  (97)  solves  Eq.  (96),  and  we  have  neglected  the  (slow) 
variation  of  5  with  We  now  consider  Eqs.  (77),  (96),  and  (97) 

in  various  limiting  regimes. 
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(a)  a  ■  a  and  A  ■  0:  For  S*y^6<<  a  ,  and  beam  radius  equal  to  the 
c  b  d  c 

conducting  wall  radius  (a  -  a^) ,  Eq.  (96)  reduces  approximately  to 
Jq(T^)  ■  0.  Making  use  of  Eqs.  (94)  and  (95),  the  real  frequency 
is  given  by  [Eq.  (85)], 

°r  -  lucz  "  t^V2  I2^  •  n*1’2 .  (98) 

r  CZ  X0ny* 

2  2 

where  Moreover,  making  use  of  ac  -  a,  B^y^S  <<  a£  and  xJq(x) 

-J^(x),  Eq.  (97)  reduces  to 

ft  1  <5  Ta  2  2 
i  *  2  ac  3ra/3nr  Vb 

(99) 

c 

and  the  slow-wave  (lower)  branch  in  Eq.  (98)  corresponds  to  instability 


(ni  >  0)  with  growth  rate 


16  „2  2i  ,1/2  “pb*®  .  „ 

0  -  —  B,  y,  VI  ,  n  »  1,2, .. . 

1  2  ac  b  b  A0nVb 


(100) 


(b)  a  »  a  and  A  ■  0:  For  6.  y,  5  <<  a  ,  large  conducting  wall  radius 

C  O  D  C 

2  2 

(a  >>  a),  and  III  a  >  1,  Eq.  (96)  reduces  approximately  to  J, (Ta)  ■  0, 

Q  1  1  f\;  X 

and  the  real  oscillation  frequency  (Jr  is  given  by  [Eq.  (83)] 


. ,,1/2  “pb*® 

(fl  -  tu  )  -  ±  u  j  ~ 

r  CZ  XlnYb 


(101) 


where  J.(X,  )  *  0,  and  use  has  been  made  of  A  ■  0  and  y  <  0.  To  the 
x  in 


accuracy  of  Eq.  (97),  we  then  find 


(102) 


for  |T|2a2  >  1. 
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On  the  oC".-;-  1: •  .i ,  for  a  »  a,  the  di sr-rs i on  re  l;«r  •  ».:i  1 

c 

i  i  2  2 

also  supports  solution.;  with  |  T ,  a“  <<  1.  In  this  cu  ~q .  (77, 
reduces  to  Eqs.  (79)  and  (81),  which  has  solutions  (for  A  =  0 
and  <<;  ac^ » 


a/2 


2£n 


(103) 


and 


.  ,  ,  ,(SV  -  ) 

„  1  J  .2  2_r _ cz 

n.  =  —  s.Yv"?  > 

1  2  ac  b  b/l  +  2  in— 

a 


\l  +  2**t)  ' 


(104) 


As  in  Eq.  (97),  the  lower  branch  in  Eq.  (102)  is  unstable  with  growth 
rate 


O  =Ii,22  I,  1 1/2  !°b 

°i  2i^bYbl^l  17 


ka 


(l  +  2*»  r) 


1/2  • 


(105) 


Comparing  Eqs.  (100)  and  (105),  it  is  clear  that  the  resistive-wall 

instability  exhibits  the  strongest  growth  when  ac  *  a. 

(c)  a  >>  a  and  A  #  0:  To  illustrate  the  stabilizing  influence  of 
c 

a  spread  a  in  canonical  angular  momentum  on  the  resistive  wall  instability, 

2  2 

we  consider  the  regime  where  | T j  a  <<  1  and  make  use  of  Eqs.  (79)  and 

(81)  to  investigate  stability  properties  for  p  <  0  and  A  i  0. 

2  2 

Expanding  Eq.  (81)  for  small  3bYb5^ac  <<:  1»  we  determine  0.^  =  Rew  and 


*  Imw  to  be  given  by 


lu 


cz 


1  5  2  2 

«i  - 1  rVb 

c 


|U!1/2.  ,ka  /  a  ' 

- - 2*-  1  +  2£n 


1/2 


jy 


2Yb 
1 1/ 2 


u  .  ka 
_ 


Yb( 1 


a^  \ 

2  Yu  U  +  j 


172 


ilki 


VRo 


(106) 


(107) 


for  the  unstable  branch  with  Imu,  ■  P  >  0.  From  Eq.  (107), 
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we  find  Chat  Che  resiscive  wall  instability  is  completely  stabilized 


by  a  small  spread  A  in  canonical  angular  momentum  satisfying 


(__* _ )2  >  1 _lL 

V^b^b^o/  16  a  2 


2  2  2  2 

>a  Vb 


|1  +  2Hn 

Of  course,  Eq.  (107)  reduces  to  Eq.  (105)  in  the  limit  A  •*  0. 


|u|c2^l  +  2in  —  j 


(108) 
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CONCLUSIONS 

In  this  paper  we  have  investigated  the  equilibrium  and  stability 
properties  of  a  relativistic  electron  ring  within  the  framework  of  the 
'linearized  Vlasov-Maxwell  equations.  The  analysis  was  carried  out  for 
perturbations  about  an  electron  ring  located  at  the  midplane  of  an  ex¬ 
ternally  applied  mirror  field  combined  with  an  applied  toroidal  field 
.  Equilibrium  and  stability  properties  were  calculated  in  Sections 
II  -  V  for  an  equilibrium  distribution  function  which  incorporates  a 
spread  in  canonical  angular  momentum  P0  [Eq.  (13)].  Equilibrium  proper¬ 
ties  were  calculated  in  Section  II,  and  one  of  the  most  important  features 
of  the  equilibrium  analysis  is  that  the  maximum  beam  density  confined  in 
the  modified  betatron  field  configuration  can  be  greatly  enhanced  by  the 

gy  J* 

toroidal  magnetic  field  Bqq  (Eqs.  (43)  and  (44)).  Moreover,  the  rota¬ 
tion  frequency  in  the  minor  cross  section  of  the  ring  plays  an  impor¬ 
tant  role  in  determining  detailed  equilibrium  properties. 

The  formal  electromagnetic  stability  analysis  was  carried  out  in 
Section  III,  and  stability  properties  were  calculated  for  eigenf requency 
u  near  harmonics  of  w  .  A  closed  dispersion  relation  (Eq.  (77)]  was  ob¬ 
tained  assuming  that  the  electron  ring  is  located  inside  a  toroidal  con¬ 
ductor  with  finite  resistivity  and  minor  radius  ac  <<  R^.  Negative-mass 
stability  properties  were  investigated  in  Section  IV  for  zero  resistivity, 
including  the  important  influence  of  equilibrium  self  fields.  For  a  low- 
density  ring,  a  modest  spread  A  in  canonical  angular  momentum  stabilizes 
the  negative-mass  instability.  In  a  gy.gh-density  ring  with  f  *  0,  however, 
if  the  self  fields  are  sufficiently  strong,  the  negative-mass  instability  can 
be  be  completely  stabilized  by  equilibrium  self-field  effects  [Eq.  (21)]. 
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The  stabilizing  influence  of  equilibrium  self  fields  on  the  negative-mass 
instability  is  an  important  new  feature  for  stable  operation  of  high-current 
modified  betatron  accelerators.  The  resistive-wall  Instability  was  exam¬ 
ined  in  Section  V.  It  was  shown  that  resistive-wall  instability  can  also 
be  stabilized  by  a  modest  spread  A  in  canonical  angular  momentum 
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APPENDIX  A 


DERIVATION  OF  APPROXIMATE  INVARIANT  P.  -  CONST.  FOR  A  CIRCULAR 


BEAM  WITH  a  ■  b  AND  n  »  1/2 


In  this  Appendix,  we  make  use  o£  Che  single-particle  equations 
of  motion  in  the  equilibrium  field  configuration  and  the  assumptions 
outlined  in  Sec.  II  to  determine  the  conditions  where 


1  eBa  a 

10  2 

PP«  *  2  “0 


(A.  1) 


is  a  good  approximate  invariant.  Here,  (p,$,0)  is  the  toroidal 
polar  coordinate  system  defined  in  Fig.  1  and  Eq.  (2). 

In  (r,6,z)  coordinates,  the  radial  and  axial  equations  of 
motion  can  be  expressed  as 


(A.  2) 


p  -  F  *  -e  (e°  +  -  v  B°  -  -  v  B°  )  . 
z  z  \z  cr0  c  0  r  / 


(A.  3) 


Here,  E®  »  -3<{>n(r,  z) /3r  and  E^  ■  -3$n(r,z)/3z  are  the  radial  and 

1  U  2  U 

axial  self  electric  field  components  determined  from  the  equilibrium 

electrostatic  potential  $Q(r,z),  and  -  BQr  +  BQr  -  -(3/3z)(AO0  +  Aq0) 

and  B®  ■  B**1  +  B®z  *  r  1 (3/3r) [r (Aq*C  +  Aqq)]  are  the  radial  and 

axial  components  of  the  total  (applied  plus  self)  equilibrium 

magnetic  field.  The  quantity  A0(r,z)  ■  Aq*C  +  Aq0  is  the  0-component 

of  the  equilibrium  vector  potential.  Making  use  of  the  assumptions 

enumerated  in  Eqs.  (3),  (4),  (5),  (7),  and  (9),  and  considering 

Q  Q/  ,2  2\ 

density  and  current  profiles  of  the  form  n/(r,z)  *  a,  +  — =•  land 

0  o/r'2  ,2\  V  t2/ 

J0^(r,z)  •  +  — 2  J*  w^*re  r'  “  r  "  B0’  *C  *s  sCra^8ht forward  to 


v  **•  «U£aJ( , ,, .  1 
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show  to  the  accuracy  required  in  the  present  analysis  that  the  equilibrium 
potentials  and  field  components  can  be  approximated  in  the  ring  interior 

<!r  -  Rni  “  lr'l  £  a  1*1  *  *>)  by8,9 


Er  “  &  r*  , 


(A.  4) 


E0  .  -Ased-f)^  ^  2  , 


(A.  5) 


+0  *  i  4ire<1*f)  ^  (br'2  +  az2)  * 


(A.  6) 


Jr  *  "nBz  R.  “  47TeR06bftb  a+b  R„ 


(A.  7) 


Bz  *  ®z  ’  nBz  IT  +  4rreRoSb&b  * 


.  (A. 8) 


rA9  *  ROA0(VO)  +  RoV’ 


1  *‘51.9 

+  ±  (l-n)B  r'“  +  ±  nB  z* 
i  Z  4  Z 


(A.  9) 


+  i  4ireRoBb  ^b (br’2  +  az2)  * 

where  r*  -  r  -  Rq,  Bz  =  »£*  (Rq.O),  ^  5  nJ(RQ,0),  Bb  s  V°(R0,0)/c, 

awf 

and  n  s  -[rHoBg^  (r,z)/Sr]£^  0)  tbe  externa^  index.  Moreover, 

from  Eq.  (1) ,  the  applied  toroidal  field  can  be  approximated  by 


B°  -  &e  (1  -  r'/R0> 


(A. 10) 


in  the  ring  Interior. 

We  make  use  of  Eqs.  (A. 2)  and  (A. 3)  and  Eqs.  (A. 4)  -  (A. 10)  to 
form  the  difference  product 
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2pr  “  r’P,  “  “7*  +  r’vr) 


+  4ire  (1  -  f)^ 


-  k  ~  a 
a  +  b 


r  z 


(A.  11) 


-  i  h.* 

c 


b  -  a  , 
i,  — 7T  r'  z 
b  a  +  b 


+  zv. 


(?4 


eB  ,  ' 
*  +  2n  —  — 

c  c  *o; 


2 

where  small  terms  of  cubic  and  higher  order  (zr'v^,  r'  vf,  etc.)  have 
been  neglected  in  Eq.  (A. 11).  Note  that  the  self-field  contributions 
in  Eq.  (A. 11)  vanish  identically  for  a  circular  beam  with  a  *  b. 


In  the  final  term  in  Eq.  (A.  11) ,  we  express  p 


6 


P0/r  +  (e/c)A0(r,z)  - 


Pn/r  +  (e/c)A„(r,z)  +  5Pfl/r,  where  <5P 


9 


0' 


0 


PQ  -  Pq,  and  make  use 


of  PQ  »  (1/2) (eBz/c)R^  and  A0(Rq,O)  »  i^Rg/2.  Neglecting  terms  of 
cubic  order  and  higher,  and  assuming  a  circular  beam  with  a  *  b,  Eq.  (A. 11) 
can  be  approximated  by 
•6. 


*p. 


*  r'P, 


8 


(zv  +  r'v  ) 
z  r 


(A. 12) 


+  zv 


fi 

9  R2  ‘ 
R0 


Strictly  speaking,  a  circular  beam  with  a  ■  b  requires  external  field 
index  n  *  1/2,  so  chat  Eq.  (A. 12)  reduces  to 


eB, 


*Pr  -  *’PZ  -  ~r  <2vz  +  r\)  +  ev0  -J- 

R0 


(A. 13) 


A-3 
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Comparing  the  two  terms  on  the  right-hand  side  of  Eq.  (A. 13),  and 

estimating  z  n,  r'  a,  vz  ^  n.  p^/y^o,  and  Ybn>ve  ^  c^Rq/c,  we  find 

that  the  6P  contribution  in  Eq.  (A. 13)  is  negligibly  small  whenever 
9 

the  inequality 


I  ".I  ^  »e  Pi 

Tb"abellO  "*  P8 


(A. 14) 


is  satisfied  [Eq.  (8)].  Within  the  context  of  Eq.  (A.14),  Eq.  (A. 13) 
can  be  approximated  by 


d_ 

dt 


d_ 

dt 


0  , 


(A. 15) 


2  2  2  2  2 

which  is  the  required  result.  Here,  p  *  r'  +  z  *  (r  -  Rq)  +  z  ,  and 
p^  is  the  4—  component  of  mechanical  momentum  in  the  toroidal  polar 
coordinate  system  (p,4,0)  illustrated  in  Fig.  1. 


A-4 


MSWC  TR  81-389 


APPENDIX  B 

ELECTRON  TRAJECTORIES  IN  THE  EQUILIBRIUM  FIELD  CONFIGURATION 

In  this  Appendix,  we  determine  the  electron  orbits  in  the 
equilibrium  configuration  described  self- consistently  by  Eq.  (13) 
in  Sec.  II  and  Appendix  A.  For  a  circular  beam  with  a  -  b  and  n  ■  1/2, 
we  express  the  Hamiltonian  H  in  Eq.  (27)  as 


H  *  +  2rs  t<p;'Ko“c0z,)2  +  (pz+Kmucer,)2] 


(6Pfl) 


(B.l) 


T5  +  "cz(1  -  r'/R0)(5V  +  7  +  z*2>  * 


where  (P',  P')  are  the  transverse  canonical  momenta  in  the  toroidal  B„ 
T  Z  u 

field. 


Pr  *  P^K^ce*'  * 


(B.2) 


p; 


2  2  2  2 

uia  *  ui  /2  +  (u>  ,/2)[9v  -  (1-f)]  is  the  betatron  frequency  defined  in 

B  CZ  pD  O 

^  A 

Eq.  (39),  and  w  *  eB  Ay.  me  and  u>  «  eB  /y.mc.  The  "primed"  orbits 

Cq  w  d  CZ  Z  0 

(r' ,z' ,8' ,p' ,p' ,p')  in  Eq.  (B.l)  pass  through  the  phase  space  point 

(r  -  Rq,  z,8,pr,pz,p0)  at  time  t'  »  t.  Moreover,  p^  -  Yb“vr  “d 

p'  *  Yvov'  for  the  nonrelativistic  transverse  motion  described  by  Eq. 

(B.l).  From  Eq.  (B.l),  Hamilton's  equations  of  motion  (d/dt')  (3H/3Pp 

-  -3H/3r ' ,  (d/dt’) (3H/3P  ')  -  -  3H/3z'  and  dB'/dt'  -  3H/3Pa,  give 

z  o 


d2r’ 


dz ' 


iP. 


V  ^75  - +  v'-cejr- +  “c,  if-  -  vv’> 


(B.3) 


B-l 
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dr1 

'b  c8  dt' 


‘  W 


(B.4) 


(B.5) 


oscillator  equations 
(B.  6) 

(B.  7) 


Introducing 

Eqs.  (B.6)  and  (B.7)  can  be  combined  to  give 


.  *  _ 


cz 


6?. 


+  iz' 


(B.  8) 


Simplifying  Eqs.  (B.3)  and  (B.4)  gives  the  coupled 

!r'  '  “c8  dtY  +  “s(t’ - . 

V  Ybm“sR0 

dr-y  +  wfz'  -  0  . 


_d 

dt 


fr6Pe)-°  • 


dt 


,2  z'  +  “ce  dt' 


“S 


d  -  v  ,  .  d  «  ,  2  f  a 

+v 


which  can  be  integrated  to  give 


n'  «  ^expCiu^t)  +  h_exp (iu^x) 


(B.  9) 


(B. 10) 


where  x  ■  t’  -  t,  ri+  are  constant  complex  amplitudes,  and  the 
frequencies  u»^  are  defined  by  Eq.  (43)], 


± 

% 


(B. 11) 


Mote  from  Eqs.  (B.8)  and  (B.10)  that  the  (r',z’)  motion  is  biharmonic 
in  w*  and  wT,  and  that  r'  and  z'  can  be  expressed  as 

D  D 


B-2 
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01. 


-•  . 


cz 


Yb"^ 


2 —  5P0  +  (r^cosu^T  +  t^cosu^T) 


-  (z^sinu^r  +  z^sinoi^T) 


(B.12) 


z'  ■  (z^cosu^t  +  S^coso^t) 

(B.13) 

+  (rJsinw^T  +  r^sinu^r)  . 


The  oscillation  amplitudes  and  z^  of  course  can  be  determined 
from  the  boundary  condition  that  the  particle  trajectories  (r' ,z' ,p^,p^) 
pass  through  (r-Rp.z.p^p^)  at  t  ■  t*  -  t  ■  0.  For  a  thin  ring  with 
a  <<  Rq,  we  note  that  |r||,  |2||  <<  Rq. 

Substituting  Eq.  (B.12)  into  Eq.  (B.5)  and  integrating  with  respect 
to  t' ,  we  obtain 


0'  -  6  +  (<ucz  -  iidPg/Y^mR^T 

-1-  (small  oscillations) , 


(B. 14) 


where  u  is  defined  by 


cz 

2 


1_ 

2 

Yb 


(B.15) 


In  Eq.  (B.14),  the  small  oscillatory  terms  (of  order  a/Rg)  are  biharmonic 


in  <d^  and  make  negligibly  small  contri'  itions  in  the  frequency  regime 

considered  in  the  stability  analysis  in  Sec.  III.  Depending  on  the 
2 

sign  of  -  (1-f)  and  the  size  of  and  “pjj/wC2»  we  note  ^rom 
Eq.  (B.1S)  that  y  can  assume  large  positive  or  negative  values. 
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